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Abstract 

The decays B® —>• and B® —>• D~ Df probe the CP-violating mixing phases 

4>d and 4> s , respectively. The theoretical uncertainty of the corresponding de¬ 
terminations is limited by contributions from penguin topologies, which can be 
included with the help of the [/-spin symmetry of the strong interaction. We 
analyse the currently available data for B® s —> Df s Df s decays and those with 
similar dynamics to constrain the involved non-perturbative parameters. Using 
further information from semileptonic B® —> Df£ + V£ decays, we perform a test 
of the factorisation approximation and take non-factorisable 5'C/(3)-breaking cor¬ 
rections into account. The branching ratios of the B® -» DfDf, B® —> D~Df 
and Bg —> D~Df, B® —> Df Df decays show an interesting pattern which can be 
accommodated through significantly enhanced exchange and penguin annihilation 
topologies. This feature is also supported by data for the B® —»• DfDf chan¬ 
nel. Moreover, there are indications of potentially enhanced penguin contributions 
in the B® —> Df Df and B® Df Df decays, which would make it mandatory to 
control these effects in the future measurements of f>d and 4 > s . We discuss scenarios 
for high-precision measurements in the era of Belle II and the LHCb upgrade. 
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1 Introduction 


CP-violating effects offer important tools to search for new physics (NP) beyond the 
Standard Model (SM). In this endeavour, B q -B q mixing (q = d, s ) is a key player. This 
phenomenon does not arise at the tree level in the SM and may induce interference effects 
between oscillation and decay processes, resulting in “mixing-induced” CP violation. The 
BaBar and Belle experiments at the e + e _ 5-factories and the LHCb experiment at the 
Large Hadron Collider (LHC) have already performed high precision measurements of 
the and B®-B° s mixing phases (ftd and (ft s , respectively [T, 21. In the era of the 

Belle II Q and LHCb upgrade |3], the experimental analysis will be pushed towards new 
frontiers of precision. 

In this paper, we present an analysis of the decays 5° —>• and 5° —>■ D~D+ 

which are related to each other through the 5-spin symmetry of strong interactions 
(5], 6 . With the help of this flavour symmetry, penguin effects can be included in the 
determination of (ftd and (ft s from the mixing-induced CP asymmetries of these decays. 
The theoretical precision is limited by non-factorisable 5-spin-breaking effects. The 
impact of these contributions can be probed in a clean way by comparing branching ratio 
measurements of the non-leptonic decays with data from semileptonic 5° —>- and 

B° s — D~£ + U£ decays. The use of the latter two modes is a new element in this strategy. 
We also explore the role of exchange and penguin annihilation topologies, which govern 
the decays 5° —>■ D~D+ and 5° —>• 0- These modes are also related to each 

other by the 5-spin symmetry of strong interactions. 

The analysis of the 5° —> D^, B® —>■ D~D+ system complements the determina¬ 
tion of 4>d and (ft s from the decays B% —> J/iftK^ and B° s —> J/iftq *), respectively, where 
penguin effects have to be included as well 15 


17 


The dynamics of the B —> DD 
modes differs from those of the B —s- J/iftX channels. I 11 the latter case, the QCD pen¬ 
guins require a colour-singlet exchange and are suppressed by the Okubo-Zweig-Iizuka 
(OZI) rule 18-20 , while this feature does not apply to the electroweak (EW) penguins, 
which are colour-allowed and hence contribute significantly to the decay amplitudes |2l|. 
On the other hand, the QCD penguins are not OZI suppressed in the 5° —> , 

5° —>• DjDf system, whereas the EW penguins contribute only in colour-suppressed 
form. The EW penguin sector offers an interesting avenue for NP to enter weak me- 
such as in models with extra Z' bosons 


son decays 22-25 


26, 27 


and would then 
should the Z' bosons have 


lead to discrepancies in the determined values of (ft d and (ft, 

CP-violating flavour-changing couplings to quarks. 

The outline of this paper is as follows: in Section [2j we discuss the decay amplitude 
structure of the 5° —>- D^D^, B ° —>- D~D+ decays and their observables, while we turn 
to the picture emerging from the current data in Section [3] There, we include additional 
decay modes, which have dynamics similar to the B® —>• 5 ° —s- D~D+ system, 

to address the importance of exchange and penguin annihilation topologies, and probe 
non-factorisable effects by means of the differential B® —>■ D^£ + U£ rate. We perform a 
global analysis of the penguin parameters a and 9 , which allows us to extract <ftd and 
c ft s from measurements of CP violation in the 5 ° — > and B Q S —> D~ Df modes, 

respectively. The current uncertainties of these measurements are unfortunately still very 
large. In Section [4| we focus on the era of the Belle II and LHCb upgrade, and explore 
the prospects by discussing different scenarios. Finally, we summarise our conclusions in 
Section |5j In an appendix, we give a summary of the various parameters and observables 
used in our analysis. 
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Figure 1: Illustration of topologies contributing to the —» D d , s) D d ,, decays. 


2 Decay Amplitudes and Observables 

2.1 Amplitude Structure 

The B d —> D d D d mode is caused by b —> ccd quark-level transitions, and in the SM re¬ 
ceives contributions from the decay topologies illustrated in Fig. [l] The decay amplitude 
takes the following form (5[: 

A ( B° d -> D- d D+) = -XA [1 - ae id e*] , ( 1 ) 


where 7 serves as a CP-violating weak phase and is the usual angle of the unitarity 
triangle (UT) of the Cabibbo-Kobayashi -Maskawa (CKM) matrix [28,29], while 


A = X 2 A[T + E + {P (c) + PA (c )} - {P w + PH W }] 


and 


ae ie = Rb 


{pH + pA(u)j _ | P (t) + PA (t)j 

T + E + {PW + PdF)} - {PW + PAW} 


( 2 ) 


( 3 ) 


are CP-conserving hadronic parameters. Here T and P ( q > denote the strong amplitudes 
of the (colour-allowed) tree and penguin topologies (with internal g-quark exchanges), 
respectively, which can be expressed in terms of hadronic matrix elements of the cor¬ 
responding low-energy effective Hamiltonian. We have also included the amplitudes 
describing exchange E and penguin annihilation PA ^ topologies, which are naively ex¬ 
pected to play a minor role 30 . However, we find that the current data imply sizeable 
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contributions for E + PAA — PA^ with respect to T + PA — P^\ The parameter 


(1 A2 ^ 

1 

V ub 

V 2 ) 

A 

v cb 


= 0.390 ±0.031 


(4) 


measures the side of the UT originating from the origin of the complex plane with the 
angle 7 between the real axis, while A = V us \ = 0.22548 ± 0.00068 is the Wolfenstcin 
parameter of the CKM matrix 31 , and A = ||/A 2 = 0.806 ± 0.017; the numerical 
values refer to the analysis of Ref. 32 . 

The [/-spin partner B° s —> D~ Df of the B® —> channel originates from the 

b —>• ccs quark-level processes. Its SM transition amplitude can be written as 


A(B«^D 7 D+)= 1 




1 ± ea'e ie V 7 


( 5 ) 


where the hadronic parameters A 1 and a!e le ’ are given by expressions which are analogous 
to those in Eqs. ([ 2 ]) and ([ 3 ]), respectively. The key difference in the structure of the 
B® —> D~Df decay amplitude with respect to Eq. ([l| is the suppression of the a'e 16 ' e* 7 
term by the tiny CKM parameter 

A 2 

e =-— = 0.0536 ± 0.0003 . (6) 

1 — X z 

Moreover, the overall factor of A is absent, thereby enhancing the decay rate with respect 
to B ° —> DjDj. Therefore, sizeable penguin effects in B% —* could be probed 

and subsequently used to estimate the penguin effects in B ° —>■ D~D+, applying [/-spin 
symmetry. 

The [/-spin symmetry of strong interactions implies the following relations between 
the hadronic parameters: 

ae ie = a!e id ' , (7) 

A = A'. (8) 

It is important to emphasise that hadronic form factors and decay constants cancel 
within factorisation in ae ie and a'e l9 , since these quantities are defined as ratios of 
hadronic amplitudes, as can be seen in Eq. ([ 3 ]). Consequently, factorisable [/-spin- 
breaking corrections to the relation in Eq. (J7]) vanish jsj 6|. On the other hand, the 
[/-spin relation in Eq. (j8j) is affected by S'[/(3)-breaking effect^] in B q —> D q form factors 
and D q decay constants (q = d,s ). We discuss these effects in more detail later. 


2.2 CP-Violating Asymmetries 

Due to B q -B q oscillations (q = d,s ), an initially present £>±meson state evolves in 
time into a linear combination of B q and B q states. CP violation in the B q —>- D q D q 
decays, which are characterised by CP-even final states, is probed through the following 

1 SU(3) symmetry refers to the symmetry group interchanging u , d and s quarks. The isospin, /7-spin 
and P-spin subgroups refer to interchanging u O d, d O s, and s it, respectively. Throughout the 
paper the mention of SU (3) refers to the /7-spin subgroup, unless specified otherwise. 
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time-dependent rate asymmetries |33|: 

I-» D-Dj )\ 2 - \A(j^t) -» -P-Opi 2 

WB,°(t) -► y-D, + )l 2 + I MB a q (t) -(• -D s -B„+)| 2 

-4cp(B,, C+)cos(AM,t) + -4gj?(g g -> g~Q+)siu(AM,t) 

cosh(Ar,t/2) + ^Ar(-S g ->■ D~D+) sinh(Ar,f/2) 


where AM 9 = — Ml 9 '* and AT, = denote the mass and decay width dif¬ 

ference between the two B q mass eigenstates, respectively. The three CP observables are 
given b\g 


A dir / p lP - n +N_ 26, ship, sin 7 
^ cp (^ ^D q D q )- , 


A^(B q ^D q D+)= Vq 
A A r(B q D~D+) = —r) q 


sin <f> q - 2 b q cos p q sin(0, + 7) + b q sin(0, + 27) 
1 — 2 b q cos p q cos 7 + 6 2 

'cos (j)q - 2 b q COS p, COs( 0 , + 7) + 6^ COs( 0 , + 27) 
1 — 2 bq cos pf cos 7 + b q 

where we have to make the following replacements for the decays at hand ( 5 |: 

je B o D - D + . 6 ■ i P s = _ ea ' e i0' . 


B S -> D- d D+ : = ne" , 


( 10 ) 

( 11 ) 

( 12 ) 


(13) 


The parameter p, denotes the CP eigenvalue of the final state and is given by +1. 
While the direct CP asymmetries Acp(B q —» D q D q ) are caused by interference between 
tree and penguin contributions, the mixing-induced CP asymmetries A^p ( B q —> D q D q ) 
originate from interference between B q -B q mixing and decay processes, and depend on 
the mixing phases 0^ and 0 S . These quantities take the general forms 

0 d = 2/? + 0* p , 0 s = -2/3 s + 0f\ (14) 


where f3 is the usual angle of the UT. The SM value of 0 S , which is given by — 2/3 s = —2A 2 p 
and hence doubly Cabibbo suppressed, can be determined with high precision from SM 
hts of the UT 


32 


= -< 2 . 092 l»S° 


(15) 


The CP-violating phases 0^ p vanish in the SM and allow us to take NP contributions 
to B° q -B°q mixing into account. 

ft is useful to introduce “effective mixing phases” 


ieff 


- r>+ 


= ^« + A K' ' 


(16) 


through the following expression |12,|16|: 

■4“p (B, D-D+) 


\fl - (^(B, -► 


= si “(7, D -D*) • 


(V) 


2 Whenever information from both B q —> f and B q —>• / decays is needed to determine an observable, 

as is the case for CP asymmetries or untagged branching ratios, we use the notation B d and B s . 
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Figure 2: Illustrations of the rescattering processes in Eqs. (23) and (24). 


D~ D~* 

where the hadronic “penguin phase shifts” A <p q q q are characterised by 


. ,D~D+ 

tan A (j) d d d = 


-2a cos 9 sin 7 + a 2 sin 27 


— 2 a cos 9 cos 7 + a 2 cos 27 

2 ea' cos 9' sin 7 + e 2 a /2 sin 27 
1 + 2 ea' cos 9' cos 7 + e 2 a /2 cos 27 
I 11 the limit a = a 1 = 0, we simply have 


= —2a cos 9 sin 7 — a 2 cos 29 sin 27 + 0(a 3 ) 


tan A (p 


D7D+ 


(18) 

= 2ea' cos 9’ sin 7 + 0(e 2 a 2 ) . (19) 


A^(B d D d Dj)\ a=0 = 0 , A^{B d -)• D d D+) | a=0 = sin j> d , (20) 

A&(B,^D-Dt)l,. o = 0, A'g?(B,^D;Dt)\a’-o = sm<t>,. (21) 

The penguin parameter ae l6 cannot be calculated reliably within QCD. Since this quan¬ 
tity is governed by the ratio of a penguin amplitude to a colour-allowed tree amplitude, it 
is plausible to expect a ~ 0.1-0.2. Applying the Bander-Silverman-Soni mechanism 34 
and the formalism developed in Refs. p35|361 yields the following estimate |6|: 


ae 


w 


| BSS 
I QCD 


0.08 x e* 205 ° . 


( 22 ) 


In the corresponding calculation, form factors and decay constants cancel as the param¬ 
eter ae lQ is actually defined as a ratio of hadronic amplitudes, as we emphasised after 
Eq. ®. However, incalculable long-distance contributions, such as processes of the kind 


B d [vr 7T ' 


,P 7T V ■ 

and 

B° -(■ {D;D+,D- d -D+, 


D~ d D* , 


/;') [K~K + , K'- K 


d;d* , 


(23) 


.] -> D- d D+ , B° [D;Dj, D’-Dt, ...]-> DJD+ , (24) 


which can be considered as long-distance penguins with up- and charm-quark exchanges 
37 , respectively, as illustrated in Fig. [ 2 J may have an impact on ae ie . 

I 11 this paper, we discuss strategies to control these effects by means of experimental 

, the penguin effects have to be taken into 


data. In the B d —> D d D d case (Eq 


account for the determination of (\) d . I 11 the B® —>■ D s D+ case (Eq. (21)), the parameter 


a' is associated with the tiny e factor and is hence doubly Cabibbo-suppressed. However, 
in view of the experimental precision in the LHCb upgrade era, also these effects have 
to be controlled. 
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2.3 Untagged Decay Rate Information 

For the analysis of the experimental data later on, it is useful to introduce another 
observable, containing the untagged rate information. It is defined as (5,6|: 


1 

A! 

2 

m Bd ^KK,%K) T Ba 

B (B d D~ d D+) 

e 

A 


m Bs ${m Dd /m Bd ,m Dd /m Bd ) r B ^ 

B(B, ^ D-Dt\ hm 


(25) 


where 

$(± y ) = v ^ 1 - (x + y ) 2 ] [1 - (x - y) 2 \ (26) 

is the well-known B —y PP phase-space function. Due to the sizeable lifetime difference 
in the I? s -meson system, y s = AT S /2T S = 0.0608 ±0.0045 1 2], a difference arises between 
the “theoretical” branching ratio defined through the untagged decay rate at time t — 0 
[5] and the “experimental” branching ratio which is extracted from the time-integrated 
untagged rate 38 . They can be related as 


39 


B{B s ^D-Dt\ heo = 


l ~y. 


.1 + ^Ar {B s —>• D s D+)y. 

= (1.09 ± 0.03) xB{B s ^ DjD+) , 


B (B s —>• D~Df ) 


(27) 


where the numerical estimate uses 


Aat(B s D~D+) = -1.41 ± 0.30 , 
extracted from the measurement of the effective B ° —» D~Df lifetime (6,40 


(28) 


The observable P[ takes the following form in terms of the penguin parameters [5]: 

(29) 


H = 


1 — 2 a cos 9 cos 7 ± a 2 


1 ± 2ea' cos 6' cos 7 ± e 2 a ' 2 


Moreover, the [/-spin relation in Eq. ([7]) implies 


H = — 


A&(B S -> D-D+) 
e L'^cp (B d —> B > d Dj) 


(30) 


Using Eq. ((7|) and keeping a and 8 as free parameters, the expression in Eq. (29) results 
in the following lower bound 41 42 


H > 


1 ± e 2 ± 2e cos 2 7 — (1 ± e) yjl — 2e ± e 2 ± 4e cos 2 7 


2 e 2 (1 — cos 2 7 ) 


7 = 73 . 2 ° 


-)• 0.908 


= sin 2 7 — 2 e sin 2 7 cos 2 7 + (D(e 2 ) 

Moreover, H allows us to put a lower bound on the penguin parameter a: 

1 + eH 


a > — 


1 -AH 


cos 7 ± 


1 + eH 
1 — e 2 p[ 


cos 7 


1 - H 
1 — e 2 H 


(31) 


(32) 


The signs have been chosen in such a way that this expression applies to the current 


experimental situation discussed in Section 3.4 
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Figure 3: Flow chart illustrating the new strategy to determine H using data from 
semileptonic B® —» D~£ + ui decays. 


2.4 Information from Semileptonic Decays 


The experimental determination of H through Eq. (25) requires information on the 


amplitude ratio |*4//k4.|, which is affected by [/-spin-breaking corrections to the relation 
in Eq. (J8]) . To avoid the limitations this brings, we propose a new method to determine 
H using data from semileptonic B q —» D~iXv^ decays, which is illustrated by the flow 
chart in Fig. [3] To this end, we introduce the ratio 


= 


T(B q -► D~D+) 


theo 


[dr(s» -> , 


(33) 


= sAVcVfiXo, 


l NF 


[1-2 bg COS p q COS 7 + bg] , 


(34) 


where the parameters b q and p q are given in Eq. (13); V cq is the relevant CKM matrix 
element, fo q denotes the D f/ -meson decay constant defined through 


(D q (p)\q^l 5 c\0) = —if Dq Pn , 


(35) 


and the factor X D is given by 


Xn = 


(m 




m 


2 \2 
D a ) 


m 2 B (m 2 R - 4 m 2 ,J 


b, 


where the form factors are defined through 


Fn qDq \ 


-1 2 


m 


D n 


F? qDq (ml) 


(36) 


(D+(k)\c^b\B"(p)) 


= Ff qDq (g 2 ) 


(p + - 


m 2 -m 2 Dq 


T 


% 


+ F 0 BqDq (q 2 ^ 


m 2 - m 2 Dq 


r 


q,, (37) 


with q = p — k. Ratios between non-leptonic decay rates and differential semileptonic 


rates as in Eq. (33) are well-known probes for testing factorisation 43 48 . 


7 
















































The parameter F measures non-factorisable effects in the amplitudes dehned through 
Eq. ([2]), which we may write as 


A g = X 2 AT, 


(9) Aq) 


fact^NF 5 


where 


71(9) • f 

d fact * Dq 




m 


B a 


m 


d c 


B n D , 


'"•'mk) 


(38) 


(39) 


is the amplitude of the colour-allowed tree topology in factorisation, with Gf denoting 
Ferrnis constant. In naive factorisation, aj^p = cq, where aq represents the appropriate 
combination of Wilson coefficient functions of the current-current operator^} We have 
A = Ad and A' = A s , and shall suppress the label q in the following discussion for 
simplicity. Introducing the abbreviations 


p(ct) ^ p(c) _ p{t) pj^ct) = PA® _ PA® 


we obtain 


where 


a NF = (1 + r P ) (1 + x) a T 


NF > 


r P = 


p{ct) 

~T~ 


(40) 

(41) 

(42) 


measures the importance of the penguin topologies with respect to the colour-allowed 
tree amplitude. On the other hand, 

PA (ct ) 


x = |x|e* CT = 


E 


(43) 


T + p( rf ) 

probes the importance of the exchange and penguin annihilation topologies. We will 


return to x in Subsection 3.3.2, where we determine |x| and the CP-conserving strong 
phase a from experimental data. The parameter a^ F describes the non-factorisable 


corrections to the “tree” diagram (Eq. (39)), i.e. we have 

T = Tfact a NF = ^fact [1 + ^Nf] 

with A^ f = 0 for exact factorisation. 

Finally, the observable H can be expressed as follows: 


(44) 


H = 

v cs 

2 

R D d 


fD s 

2 

\X D .] 


0) 

a NF 

2 

v cd 


. r d s . 


_ f D d _ 


. X D d . 


(<f) 

a NF 



(45) 


In comparison with Eq. (25), the advantage is that the theoretical precision is now only 


limited by non-factorisable 17-spin-breaking effects. Moreover, as the R Pq are ratios of 
B q rates, the dependence on the ratio of fragmentation functions f s / fd, which is needed 


for normalisation purposes 50 , drops out in this expression 


It is instructive to have a closer look at the non-factorisable 17-spin-breaking effects 
entering Eq. (45), although we will constrain them through experimental data. We 


obtain the following expression: 


(s) 

a NF 


1 4 . /\A 

1 + ^NF 


1 + r^p 


1 + X® 


1 4- A T ® 

1 + ^NF 

a {d) 

U NF 


1 A ( d ) 

1 + ^NF 


X _)_ r A> 


1 + x( d l 


1 + A T(d) 

1 ^ ^NF 


( 46 ) 


3 For a detailed discussion, see Ref. 49 . 











































Using heavy-meson chiral perturbation theory and the 1/Nq expansion, non-factorisable 
SU (3)-breaking corrections to the colour-allowed tree amplitudes of B q —* D q D p decays 


were found at the level of a few percent in Ref. 51 , suggesting small corrections from the 
last factor. The U -spin relation between the B ( ) —* D~D+ and B ° —>■ Dy decays is 
reflected by the one-to-one correspondence of the r^p and terms. These contributions 
enter Eq. (46) only in ratios of terms with structures 1 + where we expect the 


to be at most 0(0.2). Assuming S'f/(3)-breaking at the 30% level for the A^ terms 
yields a correction of only 0(5%) for the ratios, i.e. a robust situation. 

Let us now exploit experimental data to probe these effects. Using the ratio Rd 
may actually determine |a^p| with the help of the relation 


we 


,(») I — 


A NFl 


= 1 + = [l — ea' cos 6' cos 7 + O (e 2 a /2 ) ] 


R 


D s 


6ni\V cs \if D X Ds ’ 


(47) 


(s) . 

where A^p is now — by definition — a real parameter and the corrections due to the 
ea' term are at most at the level of a few percent. Assuming 


Anf - a nf t 1 - €su( 3)] 

with the SU (3)-breaking parameter £su( 3 ) oc "u/Aqcd, we obtain 

= ! + A^p £su( 3 ) ± O ^A^p 2 ) . 


A s ) 

7-NF 


7-NF 


1 + A (s) 


NF 


1 + A NF I 1 - 5sC/(3)] 


(48) 


(49) 


Consequently, the information for the semileptonic differential rate allows us to quantify 


the non-factorisable U -spin-breaking corrections to the determination of H (Eq. (45)). 

Let us now return to discuss the remaining quantities entering Eq. (45). The D q - 
rneson decay constants (q = d, s) can be extracted from leptonic decays: 


rm 


Q 2 

-»■ = gj m D q ml 


m, 


m Da 


I v cq \ 2 l 


2 f 2 

D a 


(50) 


The current experimental status has been summarised in Ref. J52|: 
f Ds = (257.5 ±4.6) MeV , f Dd = (204.6 ± 5.0) MeV , f D Jf Dd = 1.258 ± 0.038 . (51) 
A detailed overview of the status of lattice QCD calculations has been given by the 


FLAG Working Group in Ref. 53 


I 11 the infinite quark-mass limit, the following consistency relation arises 54 

nB„D, 


r 0 


(q 2 


F? qDq (q 2 ) 


= 1 


<T 


(m B +m Dq y 


(52) 


For a discussion on QCD and Aqcd /dtiq corrections to this relation we refer the reader to 


Refs. 54-56 . There has recently been impressive progress in the calculation of hadronic 


form factors within lattice QCD, where now the first unquenched calculations of the 
B —> form factors at nonzero recoil are available 57,58 . Using the results of 

Ref. 58 , we obtain 


TpB qD q 


KJ 


jpBqDq 




= 0.917 ±0.079, 


(53) 
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while the expression in Eq. (52) gives 0.924, thereby indicating small corrections. In the 
numerical analysis in this paper, we will use the result in Eq. (53). 

Experimental data for the semileptonic decay B—y Dj£ + u^ is not yet available. 
Although it is experimentally challenging to disentangle the semileptonic B® —>• D~i + v^ 
and B® —> D*~£ + U£ decays, it might be feasible to distinguish them due to the shifted 
invariant mass spectrum of the Df combinations, and the difference in the missing 
reconstructed mass, which is correlated to the “corrected mass” as illustrated in Ref. 


59 


Combined with a fit to the angular distributions, this gives information on the different 
form factors. We encourage to add this channel to the experimental agenda of the LHCb 
and Belle II experiments and perform detailed studies for the upgrade era. On the other 
hand, the differential rate of the B ° —>■ D d £ + V£ mode has already been measured and 
will actually be used in the next section to estimate the non-factorisable effects in B d 
decays. 

The lack of experimental data on semileptonic B® decays can be circumvented by 
studying the ratio of other B d and B decays, discussed in the next Section, and applying 
517(3) flavour symmetry. The SU (3)-breaking non-factorisable effects in the ratio of B d 
and B° s decays are estimated from the deviation from factorisation in B d decays. 


3 Picture Emerging from the Current Data 

3.1 Overview 

The main objective of this analysis is the determination of the B® -B° q mixing phases 
(j>d and (ft s from the B d —> D d D d and B° s —>• D~ Df channels, respectively. High pre¬ 
cision determinations of these phases require us to control not only the contributions 
from penguin topologies, but also the impact of additional decay topologies and non- 
factorisable effects. The latter two aspects cannot be quantified using information from 
the B d —y D d D d and B ° —>■ DjDf decays alone. Additional B —>■ DD decays with sim¬ 
ilar dynamics to the B d —$■ D d D d , B ° —>■ D~ D+ system therefore need to be studied. 
An overview of the different decay modes discussed in this section and their applications 
is given in Table [l] 

3.2 Preliminaries 

The direct and mixing-induced CP asymmetries of the B d —>■ D d D d decay and the H 
observable, when using the [/-spin relation in Eq. ([7]), depend on the four parameters a, 
6, (ftd and 7 . In 1999, when this decay was originally suggested by one of us, the determi¬ 
nation of the UT angle 7 was the main goal. The proposed strategy therefore assumed 
input on (pd, determined from B d —> J/'ipK g and complemented with B® —> J/'ipKg |l 6 ], 
to determine a, 9 and 7 from Acp(B d —>■ D d D d ), A™p(B d —y D d D d ) and H js]. How¬ 
ever, at present it is possible to extract 7 in a powerful way through pure B —> 
tree decays. Using current data for these channels, the CKMfitter and UTfit collabora¬ 
tions have obtained the following averages: 

7 = (73.2+^)° (CKMfitter |32|) , 7 = (68.3 ± 7.5)° (UTfit |60j) . (54) 

For the numerical analysis in this paper, we shall use the CKMfitter result. By the time of 
the Belle II and LHCb upgrade era, much more precise measurements of 7 from pure tree 
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Decay 

A 

T 

Topologies 

P E PA 

A 

Used for: 

B° d BID+ 

A 

X 

X 

X 

X 


determination of a and 9 (and (f ) d ) 

B S D' d Dt 

A 

X 

X 




non-factorisable effect a^ F 

Bj -A D-D+ 

Aepa 



X 

X 


quantify E + PA contribution x 

-A d 7 d: 

A 

X 

X 

X 

X 


physics goal (ft s 

B°, -A D;D+ 

A 

X 

X 




SU(3 ) breaking non-fact, onf/onp 

B°, -A D d D+ 

•A'ePA 



X 

X 


quantify E + PA contribution x' 

B + -> D°Dj 

A c 

X 

X 



X 

quantify A contribution ... 

B + D°Df 

A c 

X 

X 



X 

... and consistency of aNF,c/<T NFiC 


Tabic 1: Overview of the various topologies contributing to the B —>■ DD decays. The 
naming convention is indicated in the second column. 


decays will be available (see Section [I]). Using 7 as an input, we may instead determine 
c f> d and the penguin parameters from H and the CP asymmetries of B° —> | 6 |. 

The penguin parameters thus determined allow us to take their effects into account in 
the determination of (j) s from the mixing-induced CP asymmetry ^4™p x (i? s —>■ DjDf). 


3.3 Comparing B — y DD Branching Fractions 

To quantify the contributions from additional decay topologies and the impact of non- 
factorisable effects in the B® —> D^, B ° —>■ D~D+ system, we need to extend the 

decay basis to modes with dynamics similar to the B° and B ° —>■ DjDf de¬ 

cays. If we replace the spectator quarks correspondingly, we obtain the B° s —> D~D j" 
and B® —> channels. These decays are again related to each other through the 

[/-spin symmetry. However, the exchange and penguin annihilation topologies do not 
have counterparts in B® —> and B® —>■ Df, which are characterised by the 

following decay amplitudes: 


A (B» -A D-D+) = -A.4 


1 - ie'V 7 


A (BS -A D d Df) = ( 1 - y ) A! 1 + eaVV 7 


where 


A = X 2 A 
ae i§ = R 


p p( c ) _ pW 

p(u) _ p(t) 
f _|_ p(c) _ p(t) 


(55) 

(56) 

(57) 

(58) 


A 1 and a'P 6 ’ take analogous expressions. If we use the U-spin flavour symmetry, we 
obtain the following relations: 


ae 


i6 


= ~a'e ie 


A = A 


(59) 
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Annihilation (A) 



Figure 4: Illustration of the annihilation topology contributing to B + -A D°D^, s y 

Moreover, there are the charged decays B + -A D °and B + -A D°Df, which are 
again related to each other through the [/-spin symmetry. These modes also do not 
receive contributions from exchange and penguin annihilation topologies. However, there 
are additional contributions from annihilation topologies, as illustrated in Fig. [4| which 
enter with the same CKM factor as the penguin contributions with up-quark exchanges. 

3.3.1 Probing Annihilation Topologies with Charged B decays 

The decay amplitudes take the following forms: 


A ( B + -A D°D+ J = -AA 
/ \2\ 

1 + ea! r e ie ' c e i ' 1 

A c = X 2 A T c + A (c) - P® 


1 - ~a c e i6c e 


where 


A ( B + -A D°Dt) = ( 1 - y ] X 


a r e* 0c = Rh 


ptA _ pm + ^ 

t + Pi e) - P! t] 


(60) 

(61) 

(62) 

(63) 


with A c denoting the annihilation amplitude; the expressions for the primed amplitudes 
are analogous. The penguin parameter satisfies 


A0 C ~ AO 


with 


A, 


r A 


[1 + r A \ 

(64) 

_ pm ■ 

(65) 


1 1 is useful to introduce the following ratios: 

E (B.-(■£>.-£>+ B±^£>„Z)±) 


’ m Bs $(m Du /mB±,m Dd /m B ±) r B ±~ 

B(B S -A D s Dj) theo 

_m B ± <S>(m D Jm Bs ,m Dd /m Bs ) t Bs _ 

1 e(B±^B„£>J) J 


( 66 ) 
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l(B d ^D;Dt,B±^D u D}) 


m B d $(m Du /m B ±,m D Jm B ±) r B ± 

re(B d -> £>JD+) 1 

m B ± <f>{m Dd /m B ,m D Jm B ) r B 

a a a _ 

[b(B± -*■ D„Df )J 


(67) 


Using the expressions for the decay amplitudes given above yields 

~ 2 r 

_. „ ^ ^ i. A 1 — 2 a cos# cos 7 + a 2 . . 

S (B,^D;D+,B±^D U D±) = T — - J - , (68) 

A c 1 — 2a c cos t) c cos 7 + a z 

2 r ~ -1 

- 7 D n-n+ D± n n±\ ^ 1 + 2ea' cos 9' cos 7 + e 2 a /2 

^(B d ^D d Df,B±^D u Df)= -jj- - -~^ . (69) 

A c 1 + 2ea' c cos u' c cos 7 + e z a z 

If we apply the £77(3) flavour symmetry (actually the U-spin subgroup), we obtain 

Ac = A , (70) 


while the isospin symmetry of strong interactions implies 

A' = A '. 


(71) 


If we neglect the annihilation contribution in Eq. (64) and assume the same penguin 
contributions in Eq. ( 68 ), i.e. a = d c , we obtain 


KB,^D 7 D+,B±^D u Dj) 


A 


A c 


V —spin 


-+ 1 


(72) 


A deviation from unity of this ratio would therefore imply either the presence of non-zero 
annihilation contributions or large SU (3)-breaking effects through Eq. (70). In the case 
of (69), the penguin parameters are suppressed by the tiny e factor and hence play a 
negligible role. Consequently, the ratio 


l(B d ^D d Df,B±^D u Df) = 


A! 


AL 


Isospin 




(73) 


essentially relies on the strong isospin symmetry. The current experimental results com¬ 
piled by the Particle Data Group (PDG) read as follows [I]: 



B(B S D~D+) = (3.6 ± 0.8) x 10' 4 , 

(74) 


B(B d -> D d D+) = (7.2 ± 0 . 8 ) x 10" 3 . 

(75) 


B{B± -A D u D d ) = (3.8 ± 0.4) x 10” 4 , 

(76) 


B(B ± ->■ D u Df) = (9.0 ± 0.9) x 10~ 3 , 

(77) 

and correspond tcQ 

E(B S -+ D~D+, B± D u D d ) = 1.08 ± 0.27 , 

(78) 


Z(B d D d Df, B± D u Df) = 0.89 ± 0.13 . 

(79) 

4 Since B° —> D s D d is a flavour-specific final state, we simply have -4Ar(-B s —> D s D d ) 
conversion of the time-integrated branching fraction into the theoretical branching ratio. 

= 0 for the 
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For the last decay combination, we may also employ the direct measurement of the ratio 
of the relevant branching fractions (6l], which is given by 


B(B ± -> Djyj) 

B(B d -± D d Dt) 


1.22 ±0.02 ±0.07. 


(80) 


This leads to 


l(B d -± D d D+, B- -± D u Df) = 0.878 ± 0.050 , 


(81) 


which has a significantly smaller uncertainty with respect to Eq. (79) thanks to a can¬ 


cellation of uncertainties in the directly measured ratio of branching fractions. We note 
the deviation from one at the 2.4 a level, which is unexpected. 


3.3.2 Probing Exchange and Penguin Annihilation Topologies 

The current PDG results for the CP-averaged branching ratios of the B d —> D d D d and 
B® —> D~ Df decays are given as follows [l]: 


-4 


B(B d -± D d D d ) - (2.11 ± 0.18) x 10 
B(B, -» TC/J+j = (4.4 ±0.5 ) x 1(T 3 


(82) 

(83) 


In comparison with (75) and (77), the branching ratio in (83) is surprisingly small. A 
similar pattern — although not as pronounced in view of the larger uncertainties — is 
observed if we compare (82) with (74) and (76). As the B d D d D d and —* D~D+ 


decays receive contributions from exchange and penguin annihilation topologies, which 
have no counterparts in the B ° —* D~D dl B + —>■ D°D d and B d —> D d D+, B + —* D°Df 
modes, respectively (see Table [l]), it is possible that the puzzling pattern of the data is 
actually due to the presence of these exchange and penguin annihilation contributions. 

Let us first have a closer look at the ratio of the amplitudes of the B ° —y D~Df and 


B d —> D d D+ channels: 


A(B° S -» D-D+) 
A(B° -> D d Dt) 

where 


A_ 

A! 


1 + ea'e 16 ' e 11 
1 + ea'e^'e* 7 


x' = \x'\e™ = 


T' + p( rf )' 
T' + P( rf )' 

E' + PAW 
T' + PW 


+ x' 


1 + ea'e 16 ' e 11 
1 + ea'e^'e* 7 . 


(84) 


(85) 


measures, in analogy to the parameter x introduced in Eq. (43), the importance of 


the exchange and penguin annihilation topologies with respect to the dominant tree 


topology; we use abbreviations as in Eq. (40). If we neglect the terms with the penguin 


parameters, which enter with the tiny e, and introduce the SU (3)-breaking parameter 


y^'.T' + PW _ 

~T r 

T + pW/T r 

1 ^ 1 f> _)_ p(cty 

T'_ 

1 + p(cty/fr 


we obtain the relation 


A(BS -y Dl£t) , , 

MBS -± D- d D}) B + 


( 86 ) 


(87) 
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The parameter g' is only affected by SU (3)-breaking effects entering at the spectator- 
quark level. Applying factorisation, where P ^ /T' = P^ ci2> /T' (see the comment after 
Eq. (22)), we obtain 


m 


m 


_ ' ~ Bg Ds "0 V'"D S J ( ooN 

ftact _2 _2 d n . 


Pn BsDs (m 2 


■r^B,Dd 

F 0 d l 




m B s ~ m D s 

lm B m Ds 

'1 + w s (m 2 Da Y 


~£ s (w s (m 2 Ds )y 

m B d - m D d _ 

V m B d m D d 

1 + w d (m 2 Ds )_ 


Ad(w d (m 2 Ds )) _ 


(89) 


Here we have taken into account the restrictions following for the corresponding B — y D q 


form factor from the heavy-quark effective theory 46 

m B + m Dl 


rpBqDq { n 2 


(q 2 ) = 


2 y/'m B q"i Dq 


1 - 


Q 


(m B +m Dq ) 2 


£q{w q (q 2 )) 


(90) 


where £ q (w q (q 2 )) is the Isgur-Wise function with 


™ q {q 2 ) = 


m 2 ^ + m 2 Dq - q 2 
2 m B m Dq 


(91) 


Studies of the light-quark dependence of the Isgur-Wise function were performed within 
heavy-meson chiral perturbation theory, indicating an enhancement of £, s /£,d at the level 
of 5% 62 . Applying the same formalism to f B Jfo d leads to estimates for the value 
of this ratio of about 1.2 63 , which are in agreement with the experimental results in 
Eq. (51). 

Since 1992, when these calculations were pioneered, there has been a lot of progress 
in lattice QCD (for an overview of the state-of-the-art analyses, see Ref. 64]). The most 


recent result for the SU (3)-breaking effects in the form factors reads as follows 65 

^B„D a f 2 


F, 


o 


(mz 


iF 0 d d l 


m V) 


= 1.054 ±0.047 ±0.017, 


(92) 


which is in excellent agreement with the picture from heavy-meson chiral perturbation 
theory. Using this result as an input yields 


eL,= 1-078 ±0.051. 


(93) 


The error quantifies only the uncertainties related to the form factors. We cannot quan¬ 
tify the non-factorisable effects. However, as they enter only at the level of different 
spectator quarks and as already the leading SU (3)-breaking effects are small, we expect 
a minor impact. 

The ratic0 

E (B,^D-Dt,B d ^D~ d Di) 


= 0.647 ±0.049 (94) 


m Ba $(m Di /m Bd ,m Ds /m Bd ) 

B(B S — >■ D s -D+ )theo 

m Bd $(rn Da /m Bs ,m Da /m Bs ) r Bs _ 

L B(B d ^D- d Di) J 


5 In view of the large uncertainty of Eq. (28), we use Aat(B s —> D s D+) = — cos</>f M for the con¬ 
version of the untagged experimental B° s —> DjD+ branching ratio into its theoretical counterpart. 
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then takes the following form: 


'A B » -> D s D t, B d “>■ D d 7)+) = l^l 2 + 2|g'||x , | cos(a/ - d') + |± 


/12 


(95) 


thereby hxing a circle for 2 / in the complex plane. The numerical value in Eq. (94) refers 


to a direct measurement of the corresponding ratio of branching ratios 61 
For the other decay combination, we obtain 


A(B° d -t- D d D+) 


'1 - ae ie e ir 


A(B« D-D+) 

\a) 

1 — ae^e* 7 . 



with 


x = \x\e ia 


'/' _|_ p(ct) 


In analogy to Eq. (86), we introduce a parameter 

|_ p(ct) 


e=W w = 


T + PW 


+ x 

a> 

E + PA^ 

f + pw 

1 + pW/T 

l + p(ct)/f\ ’ 


ae %d e* 7 
ae^e * 7 


'T 


T_ 



(96) 


(97) 


(98) 


which is given in factorisation by 


£?fact 


- "4 

a 


m 


B, 


m 


D s 


? B d D d 


m k) 


pBs D s t 




£?fact 


= 0.928 ± 0.044 , 


(99) 


neglecting the tiny difference between the form-factor ratios for q 2 = m 2 Dd and m 2 Ds . As 
in Eq. (93), the error quantifies only the form factor uncertainties. 

The penguin parameters do not enter Eq. (96) with the tiny e. However, if we use 
the SU(3 ) relation 

p(ut) p(ut) 

( 100 ) 


T + PW T + PW ’ 

where decay constants and form factors cancel in factorisation, we get 


with 


ae id = ae iG 


r P A = 


1 + r PA 
1 + x 

PA&t) 
pM ’ 


the parameter x was introduced in Eq. (]43|). Consequently, we have 

1 


ae id e B 


1 — ae* e e* 7 


- = 1 + O(ax) + 0(ar PA ) , 


( 101 ) 


( 102 ) 


(103) 


where the second-order terms are expected to give small corrections at the few-percent 
level. Introducing the ratio 


K B d D d D d’ B s D s D a 


m Bd $(rn D Jm Bs ,m Dd /m Bs ) r Bg 


B(B d ^D- d D+) - 

m Bs <S>(m Dd /m Bd ,m Dd /m Bd ) t r ^ 


B(B, -A D-D *) am J 


= 0.59 ±0.14, (104) 
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Figure 5: Illustration of exchange and penguin annihilation topologies contributing to 
B° d -d DjDt and B« D d D+. 

we obtain 


’■( B d -> D d D d , B s £> s DJ) = |0| 2 + 2|j||i|coe(u- a) + \x\ 


(105) 


in analogy to Eq. (95). 

It is interesting to consider the double ratio 

E(B,^D;Dj,B d ^D d Dt) 

3 (B d ^D~ d D+,B,^D-D+) 


\q'\ 

I 2 + 21^1 

\x'\ 

cos (u/ — a') + | 

x'\ 

2 

\g\ 2 + 2\g\\x\ 

cos(u; — a) ± 

\x 

2 


(dfat)“ = 1-35 ±0.26, (106) 


where we have neglected the |x^'| terms and have used Eq. (99). The experimental 
results in Eqs. (94) and (|104) give 


l(B s ^DjDf,B d ^D d Dt) 

l(B d ^DjD+,B a ^D 7 D+) 


= 1.11 ±0.28, 


(107) 


which is in agreement with the expectation in Eq. (106). The current uncertainties are 
unfortunately too large to draw any further conclusions. 


3.3.3 Probing Exchange and Penguin Annihilation Topologies Directly 


The exchange and penguin annihilation topologies can be probed in a direct way by 
means of the decays B® —y D d D d and B d —y D~Df. These modes receive only contri¬ 
butions from exchange and penguin annihilation topologies (7,30 , as illustrated in Fig. [5j 
and are related to each other through the 17-spin symmetry. The current experimental 
information on the corresponding CP-averaged branching ratios is given as follows (l]: 


B[B d DjD+) < 3.6 x 10" 5 (90% C.L.) , (108) 

B(B, -> D d D+) = (2.2 ± 0.6) x 10" 1 . (109) 


The experimental signal for the B (? —> D d D d decay is in accordance with the picture 
emerging from the discussion given above. 
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Let us now have a closer look at these decays. Their amplitudes can be written as 


A (B» -> £>j£)+) = -A Aepa [1 - a EPA e‘ e ^e‘''] , 


A (B° -> D~ d D+) = ( 1 - y ) A' epa 


where 


Aepa = A 2 A 


E + PA 


1 + ea! EPA e ie ' EPA e 11 

(ct) 


appAe 


Mepa — 


= Rh 


PA 


{ut) 


, ( rf ) 


EpPA 

the primed parameters are defined in an analogous way. We obtain then 
A(B« -4 D d D+) 


A(B d A D- d Di) 


f A' epa \ 

1 ± ea' EPA e l9 ' E PAe 11 

V A' ) 

1 ± ea'e^V 6 7 


= -i'x' 


( 110 ) 

( 111 ) 

( 112 ) 

(113) 

(114) 


where we have neglected the terms proportional to the tiny e factor and introduced the 
parameter 




/' + PA 


'(ct) 


fn d m Ds 


= 0.700 ±0.042. 


(115) 


E' + PA'A) \m Dd f Ds 

In the estimate of this S'17(3)-breaking parameter, we have used that exchange and 
penguin annihilation topologies, which are genuinely of non-factorisable nature, are ex¬ 
pected to be suppressed by the smallness of the B- and H-meson wave functions at 
the origin, which behave as Jg/mB and /d/uid, respectively 66 . The fB a /tn Bs terms 
cancel in Eq. (115), and the error of the numerical value describing the leading SU( 3)- 
breaking effect corresponds only to the uncertainties of the H^-meson decay constants 
and masses. The ratic0 


’■(B,^D d D+,B d ^D d Df) 


m Ba $(m Dd /m Bd ,m Ds /m Bd ) 

B{R S B d D d ) theo 

m Bd <S>{m Dd /m Bs ,m Dd /m Bs ) r B ^ 

L B(B d ^D d Dt) J 


= 0.031 ±0.009 (116) 


takes then the simple form 

H (B, l D' d D*, B d -4 D d Df) = |c'i'| 2 , (117) 

which fixes a circle with radius \<?x'\ around the origin in the complex plane. 
Concerning the B d — > D~Df decay, we have 

v n- n+'\ / A \ Ti „ 

= SX , (118) 


A(B° S D~D + d ) 


f A E pa\ 

A - a E pAe i6EPA e il ~ 

V A ) 

i 

1 

Pi 

a 

Cb 

- 3 ' 

i_ 


where we have neglected the penguin annihilation contributions on the right-hand side, 
and have introduced 


= 


E + PA 


(ct) 


E + PAW \m Ds f D 


fo s m Da 


= 1.408 ± 0.057 ~ - . 


(119) 


6 We assume A&r{B s —> D d D d ) = —cos for the conversion of the untagged experimental 

B ° -a- D d D d branching ratio into the corresponding theoretical branching ratio. 
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Figure 6: Determination of the parameters \x'\ and a' introduced in Eq. (85) from a fit 


to Eqs. (94) and (116), which characterise the currently available experimental data. For 


the left plot, a value u' = 0 is assumed. 


As in Eq. (115), we expect that the numerical value describes the leading 577(3)- 


breaking effect (the uncertainty corresponds only to the decay constants and masses). 
Non-factorisable SU (3)-breaking contributions to this quantity cannot be estimated at 
present. 

For the comparison with the experimental data we introduce 


l(B d ^D-D+,B,^D 7 D+) 


m B d $(m D Jm Bs ,m Dd /m Bs ) r Bs 

" B(B d ^D;Di) 

m B <S>{m D Jm B m D Jm B ) t b 

L * d d d J 

_B(B S —f D s Dj) theo. 


which takes the simple form 

Z(B d -> D;Df, B, -> D 7 D+) = |ci | 2 . 

Also in this case it is interesting to consider the double ratio 


B(B d ^D 7 Dj,B,^D 7 D*) 
E (B,^D- i D*,B i ^D- i Dt) 

which allows us to test the relation 

f B«f o 


SX 


<;'x' 


X 


f B D d 


£?fact 


X . 


< 0.107 (90% C.L.) , 
( 120 ) 

( 121 ) 

( 122 ) 


(123) 


Let us now convert the experimental results in Eqs. (94) and (116) into constraints 


on \x'\ and a'. While the latter observable simply Exes a circle around the origin in the 
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Figure 7: Current experimental constraints, given in Eqs. (104) and (120), on the pa¬ 


rameters |x| and a introduced in Eq. (97). For the left plot, a value u' = 0 is assumed. 


complex plane, the former requires information about the S'71(3)-breaking parameter q'. 


We shall use the result in Eq. (93) for our numerical analysis as a guideline. A fit yields 


nt'ix) = -0.2581SS , Xm(i') = 0.0 ± 0.14 , 


(124) 


or alternatively 


\x 


_ n 9CQ+0.031 

— U.zoo _ a039 


a 


J = (180 ± 34) c 


(125) 


with the corresponding confidence-level contours shown in Fig. [6] 

The exchange and penguin annihilation topologies play hence a surprisingly promi¬ 
nent role in the decays at hand, pointing towards large long-distance strong interaction 
effects. An example of such a contribution to the exchange topology is given by 


b° s ->[d;d+,d 


d: 


->■ d-,d* , 


(126) 


as 


illustrated in the right panel of Fig. [2] 


A similar analysis can be performed for the observables in Eqs. (104) and (120), 
which allow the determination of |i| and a. In contrast to the determination of |x ; | and 
the penguin effects in the amplitude ratios do not enter with the tiny e and lead to 
additional uncertainties. In Fig. [7], we show the constraints from the current data, which 
are still pretty weak. Here we may have long-distance rescattering contributions from 
processes of the kind 


B° 


d ^{D- d D+,D‘ d -D+, 


—y d, D Q 


(127) 


In the future, following these lines, the comparison between the values of x' and x 


will offer yet another test of the relation in Eq. (123), going beyond Eq. (122) through 
information on the strong phases. 
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3.4 Global Analysis of the Penguin Parameters 

3.4.1 Information from Branching Ratios and Non-factorisable Effects 


Let us now use the currently available data to constrain the penguin parameters. Un¬ 
fortunately, a measurement of the differential semileptonic B® —> D~l + vt rate is not 


available. Consequently, we may not yet apply Eq. (45) and have to follow a different 


avenue, involving larger theoretical uncertainties. In analogy to the H observable defined 


in Eq. (25), we introduce the following quantities: 


1 

A! 

2 

m Bs $(m Dd /m Bd ,m Ds /m Bd ) 

B(B,^D 7 Dj) theo 

e 

A 


m B d &(rn D Jm Ba ,m Dd /m Ba ) r B ^ 

B [B d -3 D d Df) 


1 — 2 a cos 9 cos 7 + a 2 
1 + 2 ea' cos O' cos 7 + e 2 a' 2 


1 

A' 

2 

'$(m Du /m B ±,m Ds /m B ±y 

e 

A c 


_ <h (m Du /m B ±,m Dd /m B ±)_ 


B (B± -> D,Df) 


1 — 2 a c cos 6 C cos 7 + a; 


1 + 2 ea' c cos 9' c cos 7 + e 2 a' 2 


(128) 

(129) 

(130) 

(131) 


If we complement H with direct CP violation in B® —>■ D~ (and the suppressed CP 
asymmetry in B® —> D^Df) and H c with direct CP violation in B ± —> D U D^ (and the 
suppressed CP asymmetry in B ± —> D u Df), we may determine the penguin parameters 
(a, 9) and (a c , 9 C ), respectively. These determinations are analogous to the determination 
of (a, 9) from H and the direct CP asymmetry in P° —> (and the suppressed CP 

asymmetry in B ° — y D~D+). The hence determined parameters offer insights into the 
va and r P A parameters introduced in Eqs. (65) and (102), respectively, and allow for a 


comparison of the relative non-factorisable contributions. 

In contrast to the measurements of the CP asymmetries, the extraction of the P, H 
and H c observables from the data requires knowledge of the following amplitude ratios: 
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’ 1 + py ty /f'~ 

T’c 

A c T + P( rf ) 
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- 1 - C 

1 + P c (ct) /T c 
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1 + xA V 
1 + x _ T ’ 

(132) 

(133) 


(134) 


These quantities are governed by [/-spin-breaking effects in the ratio of the colour-allowed 
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tree contributions, which we may write as 
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( 135 ) 


(136) 

(137) 


where the parameters a^ F describe non-factorisable contributions affecting the colour- 
allowed tree amplitude (see Eq. (44)). If we assume that all the a^ F parameters are 


equal to one another due to the SU( 3) flavour symmetry, the following relation can be 
derived: 

_ 2 

rpf rji/ rjif 

(138) 

where the decay constants and form factors cancel. In terms of branching ratios, using 
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~f 
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f 
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Eqs. (132)—(134), this relation implies 


g = 


B (B ± -> D t 

pi) 



'B(B d ^D- d Di) 

B (B ± D t 

pi) \ 



B(B d ^D d D*)_ 


i. 


(139) 


The current data give 


Q = 0.85 ±0.16 , 


(140) 


which is consistent with Eq. (139) within the uncertainties. 

Using data for the semileptonic B° —y decay, the non-factorisable effects can 

be probed through 


Rd. = 


r(B° -7 D d D+) 




= 6 *r 2 | V^fiX 


2 f'2 v D a \~l 12 


B d D d ' a NFI 


1 ± 2 ea' cos 9' cos 7 ± e 2 a' 2 


where 
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and 

with 


a NF — a NF [1 + r p\ 1 


r' P = 


piety 
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(143) 

(144) 
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Figure 8: Flow chart illustrating the classic strategy to determine H using data from 
B —y DD branching ratio measurements. 


defined in analogy to Eq. (42). Experimentally, we find 

R Dd = (2.90 ± 0.41) GeV 2 , 
corresponding to a non-factorisablc contribution 


®nfI — I® 


NFI 


\l + r' P \ = 0.756 ±0.085, 


(145) 


(146) 


where Eq. (53) has been used for the ratio of form factors. In this result, the penguin 


effects suppressed by e in Eq. (141) were neglected. It is plausible to interpret the 


deviation from one at the 2.9 a level as footprints of sizeable penguin effects. We shall 
discuss this parameter below. 


The ratios in Eqs. (132)—(134) can be written in the following forms: 
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which is also graphically illustrated in Fig. 
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(147) 

(148) 

(149) 


fact 


cznf introduced in Eq. (41) through the (1 


In the above equation, differs from 


x ) term, to account for the contributions 
from exchange and penguin annihilation topologies, which are absent in the decays of 
the other two ratios. We thus have 


®nf — (1 ± x)a 


(0) 
NF ’ 


(150) 
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As in the discussion in Subsection |2.4 it is convenient to write 

|onf| = 1 + A' nf [1 — £su(3)\ , 


I a NF I 


= i + A; 


NF > 


( 151 ) 


where the parameter £,su( 3 ) describes SU (3)-breaking effects in the non-factorisable con¬ 
tributions, yielding 


NF 


a N F 


1 + A 7 nf 
1 +A 


NF 


— 1 + €su( 3 )A' nf + C>(A' nf ) . 


(152) 


Consequently, A' NF = —0.244 ± 0.085, as determined from Eq. (146), allows us to take 


non-factorisable corrections to Eqs. (14T)—(149) into account, assuming 



a (0) ' 

U NF 


a NF 


a NF,c 



(0) 

a NF 


Unf 


<5nf,c 



(153) 


In these relations, SU( 3)'-breaking effects enter only through different spectator quarks 
and are expected to be small. 

In order to determine the SU (3)-breaking effects in the ratio of the colour-allowed 


tree amplitudes |T'/T|, we use again Eq. (90) to derive the following expression |5]: 

= 1.356 ± 0.076 . (154) 
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For the calculation of the numerical value, we have used Eq. (92) and the values of the 


decay constants in Eq. (51). In analogy, we get 
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= 1.167 ±0.066. (155) 


In the case of the charged decays, a particularly simple situation arises in factorisation 
as the form factors cancel: 

T' f_D. 

fo d ’ 


T r 


(156) 


fact 


as is also evident from Eq. (137). 


Using the amplitude ratio in Eq. (147), we can now determine the observable H (see 


Eq. (25)). For this, we also have to quantify the uncertainty from the SU (3)-breaking 


corrections to the term involving the x and x' parameters. The analysis discussed in 


Subsection 3.3.2 allows us to accomplish this task. Using the results for x' and a' in 


Eq. (125) and the relations 


/ x 
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m B s m Ds 
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yields 


1 + x' 


1 ± X 


= 0.930 ±0.020 


(157) 


(158) 
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Figure 9: Pie charts illustrating the uncertainty budget of the H observables. 
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Finally, we obtain 


and correspondingly 


A' 

A 

A' 

A 

s 

A c 


1.261 ±0.091, 

1.167 ±0.081, 

1.258 ±0.063, 


(159) 

(160) 

(161) 


H — 1.30 ± 0.26 , (162) 

H = 1.28 ±0.29, (163) 

H c = 1.18 ± 0 . 21 . (164) 


In Fig. [9] we show the corresponding uncertainty budgets. Should the relations in 
Eq. (153) actually receive corrections, they would affect this error budget. I 11 the future, 


implementing the strategy proposed in Section 2.4, this approximation/assumption is no 
longer needed. 

The result on H allows us to put first constraints on the penguin parameters a with 


the help of the lower bound in Eq. (32): 


a > 0.052 , 


(165) 


where we have used the lower value of H at one standard deviation. 


3.4.2 Information from CP Asymmetries 


Let us now add experimental information on CP violation to our analysis. Concerning 
the decay B d —» D d D d , the current status of the measurement of the direct and mixing- 
induced CP asymmetries is given as follows: 


^CP (B d -A D- d D + d ) 


-0.07 ± 0.23 ± 0.03 (BaBar |67|) 
-0.43 ± 0.16 ± 0.05 (Belle | 68 j) , 


(166) 


A^(B d D d Dj) 


±0.63 ± 0.36 ± 0.05 (BaBar |67|) 
±1.06loii ± 0.08 (Belle | 68 :j) . 


(167) 


The measurements by the BaBar and Belle collaborations are not in good agreement 
with one another, in particular for the mixing-induced CP asymmetry. HFAG gives the 
following averages | 2 |: 


-► D- d D+) = -0.31 ± 0.14 , A%?(B d ^D- d D+) = 0,98 ±0,17, (168) 


which have to be taken with great care 
as input for the strategy discussed above 

= 0.028 for 4 


Eq. (|162j) yields x 
solution 


2 
min 


TZe[a] = 


-0 29 +0,27 

u - zy —0.20 i 


It is nevertheless interesting to use these results 
A x 2 fit to Eq. (168) and the value of H in 
degrees of freedom (a,9,(j) d , 7 ), and results in the 

Im[o] = -0.204 ;°;?o5 , (169) 
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corresponding to 

a = 0.351S , 0 = (215«‘)° , 


(170) 


and 



result in the penguin phase shift 


A0 ( 


D ~ d D l 

d 


= (sol?-'' 


32 / 


(172) 


For the fit, we use the expressions for the CP asymmetries in Eqs. (10) and 0 , and the 
expression for H in Eq. (29) with the [/-spin relation in Eq. (J7|). Since a' enters Eq. (29) 
with the tiny e, [/-spin-breaking corrections to Eq. (JtJ) have a very minor impact. 

Although the uncertainties are large, the results from the fit may indicate significant 
penguin contributions. Should this actually be the case, long-distance effects, such as 
those illustrated in Fig.pl would be at work in the penguin sector of the Bjj —>• D~D+ de¬ 
cays. It is interesting to have a closer look at the result of non-factorisable contributions 


to B d —y D d D+ in Eq. (146). It corresponds to 


|1 + f ' P | = (0.750 ± 0.055)/|5nfI 


(173) 


Assuming that |a^ F | = 1, i.e. that the colour-allowed tree topologies have negligible 
non-factorisable contributions, and that all parameters are real, this results in 


f'p = -0.250 ± 0.055 


(174) 


Note that the uncertainty on r' p only reflects the uncertainty on the input quantity 


(146), but does not take into account further theoretical uncertainties associated with the 


made approximation. Further assuming p( ut Y = p( ct ) , ; which leads to the approximation 
ae lB ~ Rbf P /( 1 + r'p), we obtain 


a = 0.130 ±0.039, 6 ~ 180 c 


(175) 


which is in the ballpark of the theoretical estimate in Eq. (22). Consequently, it is still 


premature to draw definite conclusions on anomalously enhanced penguin contributions 
at this point and future analyses are required to shed light on this issue. 

Large penguins would have important consequences for the determination of the B° s - 
B° s mixing phase (j) s from the mixing-induced CP violation of the B° s —y D~Df decay. 
The LHCb collaboration has recently presented the first analysis of CP violation in this 


channel 69 , yielding 


-eff 
s,D~D2 


3rD + = (l-l ±9-7 ±1.1)° , 
which can be converted into 


A +1 = 0.91 ±0.18 ±0.02, 


(176) 


A&(B s ^D~D+) = 0.09 ±0.20, 


^ X (5 S ->• D~D+) = 0.02 ± 0.17 . (177) 
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Figure 10: Illustration of the determination of the penguin parameters a and 6 from a 
X 2 fit to the CP asymmetries of the decay B ° —>■ and the observable H. 


These results are in good agreement with the predictions based on a global SU (3) analysis 


of the B —>■ DD system 


70 


Moreover, we obtain 


Aaf{B s D~D+) = -0.995 ± 0.019 , 


(178) 


which should be compared with Eq. 
effective lifetime +. 

U s U s 


corresponding to the direct measurement of the 


as 


If we generalise the [/-spin relation in (JT|) 

a' = £a , 9' = 6 + 8 , 


(179) 


with the assumed [/-spin-breaking parameters £ = 1.00 ± 0.20 and 8 — (0 ± 20)° (which 
are of similar size as the corresponding parameters in Ref. fl6j) and use the expression in 


Eq. (19), the penguin parameters in Eq. (170) determined from the fit can be converted 
into 


A 4>° 7Dt = - (1.73:1 (stat)«;? ((/-spin))' 


(180) 


Finally, we can extract (f> s form the effective mixing phase in Eq. (176), which yields 

4> s = ~ (0.6+gJ (stat)l °:7 (U- spin))° . (181) 

Despite the suppression through the parameter e, penguins may have a significant impact 
on the extraction of 4> s and have to be taken into account. This will be particularly rele¬ 
vant for the LHCb upgrade era. In this new round of precision, we will also get valuable 


insights into the validity of the [/-spin symmetry, parameterised through Eq. (179). 

Unfortunately, there is no measurement of CP violation in B J? —y D~ available, 
which would be very interesting, in particular in view of the situation for B® —y . 

Consequently, we may not yet determine a and 6 from the data. 

However, for the charged B + —y D decay, the PDG gives 


AcpiB* -> D u Dj) = + 0.03 ± 0.07 . 


( 182 ) 
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Figure 11: Illustration of the determination of the penguin parameters a c and 9 C from a 
X 2 fit to the direct CP asymmetries of the decay B + —> D D J" and the observable H c . 


If we complement this measurement with the value of H c in Eq. (164), we may perform 


a fit to the charged decays. It has three degrees of freedom (a, 0, 7 ), and results in the 
solution 

Tie [a c ] = —0.22lo(i8 , Xm[a c ] = 0.018 ± 0.044 , (183) 

which corresponds to 


a c = 0.22 


+0.18 
- 0.22 > 


e c = (175 ± 12) c 


(184) 


In Fig. 11, we illustrate the corresponding situation, which complements Fig. 10 

It is interesting to compare the penguin parameter ae l9 with its charged decay coun¬ 
terpart a c e 9c . We obtain 


= 

' PA — 


1 


a c e 9c 

1 + r A 

1 + X 


ae ie 

1 + r PA 


(185) 


where we have used Eqs. (64) and (101). The precision that can be obtained with the 
current data does not yet allow us to draw any conclusions regarding rp A . However, in 
the future it will be interesting to monitor this quantity as the experimental precision 
improves. Moreover, a measurement of the direct CP violation in the B ° —» D~D j" 
channel will allow us to determine ae l 9 from the information from H. The comparison 
with a c e 9c will yield the r A parameter from Eq. (64), so that Eq. (185) will then allow 


the determination of rp A ■ Consequently, following these lines, we may reveal the impact 
of the annihilation and penguin annihilation topologies in the decays at hand. 


4 Prospects for the LHCb Upgrade and Belle II Era 

Let us conclude the discussion on the B —> DD decays by exploring the potential of 
these decay modes in the Belle II era and at the LHCb upgrade. We do this using several 
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scenarios, examined in Section 4.3, that reflect the different possibilities still allowed by 


the current data. The inputs used in these scenarios are discussed first. Section [4T] gives 
the experimental prospects for the r elevant CP and branching ratio information of the 
B —y DD decays, while Section 4.2 deals with the future constraints on the additional 
decay topologies. 


4.1 Extrapolating from Current Results 

The 5-factories have pioneered the study of B —> DD decays, including the discoveries of 
numerous B —>■ DD decay modes 


71 


the measurements of branching fractions 


and the analyses of CP asymmetries 67,68,72 . The LHCb collaboration subsequently 


72 


73 


continued the study of B —* DD decays, notably focusing on the analysis of B° s de- 
which are abundantly produced at the LHC. Based on the successful 


40,61,69 


cays 

performance of LHCb during run I of the LHC, an estimate can be made of its per¬ 
formance with the data samples that are expected to be collected after the upgrade of 
the LHCb detector. For these extrapolations, an integrated luminosity of 5 fb _1 in run 
II, from 2015 until 2018, is assumed. In addition, the B production cross section will 
increase at a centre-of-mass energy of 13 TeV compared to 8 TeV by about 60%. For the 
upgrade scenario, an integrated luminosity of 50 fb _1 is assumed with increased trigger 
efficiency, leading to about a three times larger data sample per fb -1 compared to the 
B yield per fb 1 at run I. Similarly, a prognosis can be made for measurements at Belle 
II, which is expected to start taking data in 2018. Here we assume that 50 times more 
data will be collected than currently is available (1 ab _1 ) . 

The expectations for the CP asymmetry parameters are listed in Table [2j The extrap¬ 
olations are done for the currently available measurements only; no attempt is made to 
forecast the precision on yet-to-be-performed analyses. For example, the LHCb collabo¬ 
ration will also determine the CP asymmetries of the 5° —)■ D^ decay, but it remains 
to be seen what the accuracy will be in comparison with possible Belle II results. 

The expectations for the branching fractions are listed in Table [3] The “current 
measurement” column reflects the best available knowledge at this moment, which in 
some cases could have been more precise if the ratio of branching fractions were deter¬ 
mined directly, rather than dividing the individually measured branching fractions. In 
the extrapolations it is assumed that the ratios of branching fractions are determined. 
Moreover, it is assumed that the systematic uncertainties due to f s /fd (4.7%), due to 
the 5-meson branching fractions (3.9%, 2.1% and 1.3% for Df, and D° mesons, 
respectively) and due to the different 5° lifetimes (2.9% (1.5%) for a CP (flavour) eigen¬ 
state), remain the same. We assume that the total experimental systematic uncertainty 


Observable 

Current Measurement Upgrade Experiment 

Af P (B d D- d Di) 
•4SK B d -> D d D+) 
Af P (B± -> £>„£>*,) 
<ff (B, -> D 7 D+) 

-0.43 ±0.16 ±0.05 68 

1.06 ± 0.08 68 
0.00 ±0.08 ±0.02 72" 
(1±10±1)° 69 


±0.05 

±0.08 

±0.02 

±2° 

Belle 

Belle 

Belle 

LHCb 


Table 2: Experimental prospects for the currently available CP asymmetry measurements 
of B —> DD decays from Belle and LHCb. 
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Obs 

Decay Ratio 

Current Measurement 


LHCb Uncertainty 



of ratio of BRs 

2011 

Upgrade 

H 

B" D- d Di/B<i -» D-Dt 

0.048 ± 0.007 

1 


14%(12%) 

8% 

H 

B°, D;D+/B % -» D d Df 

0.050 ± 0.008 ± 0.004 (61 


18% 

7% 

H c 

B+ ->• D°Dj/B + D°D+ 

0.042 ± 0.006 

1 


15% (7%) 

6% 

z; 

B°, D-Dt/Bl -» D~ d Dt 

0.56 ±0.03 ±0.04 | 

T 


9% 

7% 


B° d -> D d D+/B« -> D-D* 

0.59 ±0.14 

l 


24% (20%) 

6% 


B° -> D- d D*/B« -> D d Dt 

0.031 ±0.009 

l 


24% (20%) 

11% 


D;D}/B» D;D* 

Not observec 





Tabic 3: Experimental prospects for ratios of branching fractions. The second and 
fourth ratios are obtained from direct determinations of the ratios of branching fractions, 
whereas the others are calculated from individual branching fractions. The value in 
brackets indicates the possible uncertainty if this ratio were determined directly. Note 
that for the calculation of the H observables, additional uncertainties due to (M//MI 
arise. 


Observable Current Measurement 

Upgrade 

fD d 

fn s 

7 

(204.6 ±5.7 ±2.0) MeV 
(255.5 ± 4.2 ± 5.1) MeV 
(68.3 ±7.5)° 

74 

76 

60 


±3.0 MeV 
±3.6 Me\ 
±0.9° (3 

75 

* 


Table 4: Experimental prospects for external input. A * indicates an average of the 
extrapolated measurement with the current PDG average. 


will decrease from 5.0% to 4.0%. In some ratios, the uncertainty on the D branching 
fractions cancels with their contribution to the f s /fd uncertainty. This is taken into 
account where appropriate. In our upgrade era scenario, systematic uncertainties will 
be the limiting factor on the ratio of branching fractions. Therefore, we would like to 
encourage research into f s /fd , —> K + K~tt + )/B(D^ —* A' - 7r + 7r + ) and B lifetime 

differences. If these three factors could be reduced to a level of about 2%, then that 
would lead to a systematic uncertainty of (5-6)% for all of these decays, assuming an 
experimental uncertainty of 4%. Finally, the prospects for improvements on external 
input parameters are listed in Table [4j 


4.2 Exchange and Penguin Annihilation Contributions 

For the construction of the H observable based on the ratio of hadronic amplitudes in 
Eq. (147), the contributions from exchange and penguin annihilation topologies, rep¬ 
resented by the parameters x' and d', need to be quantified. With future, improved 
measurements of the B —>■ DD branching ratios it is expected that the picture emerging 
from the current data, discussed in Section 3.3.2, can be sharpened further. We therefore 
explore the precision that can be achieved towards the end of the Belle II era and of 
the LHCb upgrade. Based on the best fit solution obtained from the current data, i.e. 
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Figure 12: Illustration of the determination of exchange and penguin-annihilation con¬ 
tributions through the parameter x (Eq. (43)) in the Belle II/LHCb upgrade era. The 
precision should be compared to Fig. |6j 


Eq. (125), and the prospects in Table [3j we start from the following input measurements: 



5(H S ->■ D s D +, B d —>• 

D~D+) = 0.673 ± 0.043 , 

(186) 


-> D d Dj,B d ->■ 

D d D+) = 0.033 ± 0.004 . 

(187) 

A fit yields 

Re[x'] = -0.258“ , 

Im[x] = 0.00 ±0.15 , 

(188) 

corresponding to 

x' = 0.258“ , 

a' — u' — (180 ± 15)° . 

(189) 


The associated confidence-level contours are shown in Fig. 12. Note that at first sight 


these uncertainties do not seem to have improved significantly with respect to the present 
experimental situation. However, this is merely caused by the shape of the confidence 


contour in Fig. 12 


4.3 Future Scenarios 


To achieve the smallest theoretical uncertainty on the H observable, it should be con¬ 


structed from the semileptonic decay information, see Eq. (45), as explained in Section 
2.4| This method is preferred over the direct ratio of hadronic branching fractions in 


Eq. (25), as it does not rely on form factor information, and is not experimentally limited 
by fs/fd ■ However, as the necessary information on dr/dg 2 (R(? —> D~l + Vi) is currently 
not yet available, we also do not have any estimates for the precision that can be achieved 
at LHCb or Belle II. For the following discussion, we will thus, like for the fits to the 
current data, rely on the original definition using the ratio of hadronic branching frac¬ 


tions, Eq. (25). 
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Using the currently available data, we have illustrated in Section 3A that it is possible 
to simultaneously determine the penguin contributions and the B d -B d mixing phase 0^ 
using the CP and branching ratio information of the B d —y D d D d decay. However, as 
the result of 0^ in Eq. (171) shows, the precision on 0^ is very limited. Also for the Belle 


II and LHCb upgrade era it will be challenging to reach precisions below the (10-20)° 
level. A high-precision determination of 0^ using the B d —)■ D d D d decay is only possible 
if the direct CP asymmetry and the H observable together are sufficient to unambigu¬ 
ously pin down the penguin parameters a and 6. In such a situation, the phase 0^ can 
then be determined from the mixing-induced CP asymmetry. As Fig. [10] illustrates, a 


and 6 cannot precisely be determined. The situation would arise either for very large 
values of a, which looks unrealistic, or if the H observable can be determined with a 
precision of well below the 5% level, which in view of the prospects in Table [3] is unlikely. 
The determination of 0^ from the B d —> D d D d decay is therefore not competitive with 
the results from the B d —> decay. It is thus advantageous to use external input on 

0^ for the analysis of the B —>■ DD decays. This extra information breaks the ambiguity 
that is still present in the confidence-level contours shown in Fig. 10, and can therefore 
improve the precision on a and 6, and thus also on A 0 S . For the future benchmark 
scenarios we therefore only focus on the high precision determination of 0 S . 


Using external input for <pd in principle makes one of the three observables used 
by the fit (^4cm or H) superfluous. As the H observable receives corrections 

from possible [/-spin-breaking effects through |„4 / /.A|, which result in large theoretical 
uncertainties, it is preferred to determine the penguin parameters using information 
on the CP asymmetries only, omitting H. Such a determination is theoretically clean. 
This will ultimately lead to the highest precision on a, 6 and A0 S . In this situation, the 
branching ratio information can instead be used to gain insight into the hadronic physics 
of the B d —> D d D d , B—> D~Df system. We can then follow the opposite path where 
the fit results for a and 6 can be used to determine the H observable with the help of 


Eq. (29), labelled 77( a ,e) below. Since a' enters there with the tiny e, the U-spin-breaking 
corrections affecting the [/-spin relation in Eq. (179) have a very minor impact. Since 


we now know the value of H, the relation (25) can be inverted to instead determine the 
ratio of hadronic amplitudes 


A! 


A 


\ 


eH t 


(a, 9) 


m B $ (m D - /m B m D +/m B ) t b o 

- LJ s a da d d 

m Df /m B ) r B ; 


B(B,^ D;Dt) theo 


B [B d -» D- d D* d ) 


(190) 


from the measured ratio of branching ratios. This experimental measurement of the ratio 


of hadronic amplitudes can be compared with the theoretical result in Eq. (159). This 


favourable strategy is illustrated by the flow chart in Fig. 13 


However, the ideal scenario described above cannot always be realised. When the 
value of the mixing-induced CP asymmetry is compatible with 1 (at the la level), its 
power to constrain the penguin parameters a and 6 is limited. This can best be il- 

In this situation, the 


lustrated using the contour plots, like Fig. [TO] or Fig. [17] below, 
annular constraint originating from the mixing-induced CP asymmetry becomes a closed 
disk, leading to a large overlap region with the direct CP asymmetry constraint. Conse¬ 
quently, it is not possible to conclusively pin down a and 6 in such a situation. Additional 
information is thus needed to improve the picture, and reach our target of matching the 
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i i Quantity is Theoretically Clean i i Affected by C-Spin-Breaking Effects 


Figure 13: Flow chart illustrating the favourable strategy to control A (j) s , which only 
requires information on the B d —> D d D d CP asymmetries. 


foreseen experimental precision on 0 S with an equally precise determination of A0 S . In 
this situation, the H observable forms an essential ingredient in the fit, and it can there¬ 
fore not be used to experimentally constrain the ratio of hadronic amplitudes. This less 


favourable strategy is illustrated by the flow chart in Fig. 14 


Given the current experimental situation, either of the two situations sketched above 
can still be realised, depending on the future world average for A^{B d —> D d D d ). To 
demonstrate the variety of situations in which we may ultimately end up, we made an 
overview of six different scenarios, covering both situations. Scenarios 1-3 represent the 
favourable situation in which the H observable can be omitted from the fit, in which we 
can determine a and 9 in a theoretical clean way, and get experimental access to the ratio 
of hadronic amplitudes. Scenarios 4-6, on the other hand, fall in the second category 
and do require information on H to conclusively pin down a and 6. All six scenarios 
are chosen to be compatible with the current experimental situation, with scenario 5 
representing the current best fit point, and have a < Rb , which is suggested by Eq. ([3]). 
Although it is mathematically possible for a to be larger than Rb, it would imply that 
the penguin topologies are larger than the tree contribution, which seems very unlikely. 
The condition a = Rb thus serves as a naturally upper limit for the size of the penguin 
contributions. The different scenarios we consider, and the resulting input values of the 
three observables (.A^p, *4 cp an d H) are listed in Table [5j The choice of input points 
can be compared with the current fit solution for a and 9 in Fig. 15, and with the current 


measurements of the B d —>■ D d D d CP asymmetry parameters in Fig. 16 


For each of the six scenarios, the individual constraints coming from Aq 
H are illustrated in Fig. [TTl The three left-most plots represent the favourable situation, 


c £, A^ x and 
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i i Quantity is Theoretically Clean | i Affected by 17-Spin-Breaking Effects 


Figure 14: Flow chart illustrating the less favourable strategy to control A cj) s , which re¬ 
quires information on the H observable in addition to the B$ —>■ CP asymmetries. 


No. 

a 

6 [deg] 

A dir 
M. cp 

Jmix 

*H C P 

H 

1 

0.15 

260 

-0.27 

0.70 

1.04 

2 

0.10 

200 

-0.06 

0.80 

1.07 

3 

0.25 

320 

-0.32 

0.35 

0.95 

4 

0.20 

230 

-0.26 

0.82 

1.12 

5 

0.34 

216 

-0.30 

0.91 

1.29 

6 

0.35 

190 

-0.09 

0.97 

1.34 


Table 5: Penguin parameters and observables corresponding to the six different scenarios. 
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Figure 15: Distribution of six scenarios in the TZe[a]-Xm[a] plane [Left] and in the 9—a 
plane [Right]. Superimposed are the confidence-level contours from the fit to the current 
data. The red circle and horizontal line represent the natural upper limit a = Ri. 


while the three right-most plots fall in the second category. For the three right-most 
plots the constraint from the mixing-induced CP asymmetry is more disk-like as the 
central value of ^4™p x is closer to one. As a consequence, the overlap with the direct CP 
asymmetry, which forms a narrow band in all cases, is too large to pin down a and 9. In 
this respect, scenario 4 should be seen as a limiting case; information of H is not strictly 
necessary if one is only interested in the 1 sigma results. 

For each of the six scenarios we also performed a y 2 fit, similar to the one described 
in Section T4, but including cj) d as a Gaussian constraint. The fit results for a and 9, 
and the associated values for the shifts Acp d and A(p s are listed in Table |6j [/-spin- 
breaking effects, parametrised by Eq. (179), have been included in the results for A0 S . 
The associated confidence-level contours are shown in Fig. [TT| In all cases we succeed 
in our goal of matching the foreseen experimental precision on </> s , see Table [2j with an 
equally precise determination of A(p s . This is also the case for and A(j) d , as illustrated 
in Table [TJ For the first three scenarios, which do not include the H observable in the 
fit, the resulting solution for 77( a ,e) and the values for the ratio of hadronic amplitudes 
are listed in Table [8] The resulting uncertainties are about a factor two smaller that 
the current theoretical uncertainties derived within the factorisation framework, and 
of comparable size to the experimental precision that can be obtained on the ratio 
of hadronic amplitudes describing the B® —> J/^K^ and B° s —> J/'ijjK g decays 16 


Consequently, the experimental determination of | ^4.'/is yet another interesting topic 
for Belle II and the LHCb upgrade. It will provide valuable insights into possible non- 
factorisable [/-spin-breaking effects and the hadronisation dynamics of the B —>■ DD 
decays. 
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1.4 



Figure 16: Distribution of the six scenarios in the pl an ± which can be com¬ 

pared to the current BaBar [Green] and Belle [Blue] measurement and the world average 
[Red], 


No. 

H 

a 

9 [deg 


A (j) d [deg 

A0 S [deg] 

1 

No 

o i 50+°- 032 

u.iou_ 0 Q 29 

260.0±1®;1 

3-6 tH 

- 0.15l£^(stat)lg;jjjj(tf-spin) 

2 

No 

n inn +0 079 

u.iuu_ 0 063 

200.01™ 

io.i±?:g 

-•0.55lg;^(stat)lg;^(£7-spin) 

3 

No 

o 2 c sn +0 036 

u.zju_ 0 037 

320.0±?;1 

-21.6l{j;? 

1 - 12 ^o:26( s tat)lg;S(£7-spin) 

4 

Yes 

0 200 + °° 63 
U.ZUU_ 0 .052 

230.01?^ 

14.4±f;| 

—0.75lo;f 3 (stat) IqJ 9 ([/-spin) 

5 

Yes 

0.340«;™ 

2i6.0+^;° 

9Q 1 + 14 - 7 
zy.±_12 i 

^ 1 .62«;^(stat)l°;l(C/-sp in ) 

6 

Yes 

o.3501o;ii5 

190.0lgi 

33 6 +15 ' 8 

OO.U_ 14 0 

—2.041^69 (stat) 1 q!32 ([/-spin) 


Table 6: Fit results for the different scenarios. 


No. (j)f [deg] A (p d [deg 


1 

46.8 ±6.9 

3.6±J;1 

2 

53.3 ± 7.5 

10.1±f;g 

3 

21.6 ±5.1 

~21.6l{jj 

4 

57.7 ±8.8 

UAt 9 7 8 7 

5 

72.3 ± 15.7 

9Q 1+ 14 - 7 
za.i_i2 i 

6 

76.8 ± 19.6 

oo c+15.8 
OO.U_ 14 o 


Table 7: Comparison between the effective mixing phase and the shift A (fid- 


No. 

H( a ,e) 

\A!/A\ 

1 

1.038 ± 0.039(a, 6) ± 0.002(£, 5) 

1.163 ±0.048 

2 

1.067 ± 0.053(a, 6) ± 0.001(f, S) 

1.179 ±0.053 

3 

0.946 ± 0.014(a, 9) ± 0.002(£, 5) 

1.110 ±0.042 


Table 8: Constraints on the ratio of hadronic amplitudes for those scenarios where the 
H observable is not needed. 
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Scenarios 1-3: 


Scenarios 4-6: 


a = 0.1S01SX 



i—i -^CP {B d -> D d D d ) 


a = 0.200«S 

e = ( 260 . 031 )" 



i—i -^cp D d D d ) 

0.8 

e = (230.0™)" 

M? b: = -(0.1518$!)" 



o H(B, -> D-D*) 

A(4f rfl; = -(0.75«;1)" 

E2 39 % C.L. 




E521 39 % C.L. 

□ 68 % C.L. 



0.6 

IZZI 68 % C.L. 

L-r: 90% C.L. 

\k. 



cr: 90% C.L. 


•^cp D d D d 

A™(B d -> D d D+) 
H(B q D~D+) 




TZq[o\ 7Ze[a\ 

Figure 17: Illustration of the determination of the penguin parameters a and 6 for the 
scenarios introduced above. The three left scenarios have favourable values for Aq p- 
^4™“, while scenarios 5 and 6 at the right bottom would require knowledge from H. 
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5 Conclusions 


In this paper, we have presented a detailed study of the system of B —> DD decays, 
exploring the picture emerging from the currently available data and proposing new 
strategies for the era of Belle II and the LHCb upgrade. We find that patterns in the 
current branching ratio measurements can be accommodated through sizeable contribu¬ 
tions from exchange and penguin annihilation topologies, which play a more prominent 
role than naively expected. This feature suggests that long-distance effects of strong 
interactions are at work, which cannot be understood within perturbation theory. 

Using data for the differential semileptonic B® —> D^£ + ug rate, we have determined 
non-factorisable contributions to the —> D^Df decay, finding effects at the 25% level, 
including certain penguin contributions. We have pointed out that a future measure¬ 
ment of the semileptonic B ° —> Dj£ + U£ channel would allow an optimal determination 
of the observable H from the B ° —>■ and B ° —> DjDf branching ratios, with 

uncertainties given only by non-factorisable [/-spin-breaking corrections, which can be 
further quantified through the comparison of the differential B j? —> D~£ + P£ rate with the 
Bg —> DjDf rate. Experimental studies of the B Q S —> D~£ + V£ modes are encouraged. 

In view of our new insights into the prominent role of the exchange and penguin anni¬ 
hilation topologies, the penguin contributions may also be more important than naively 
expected. The current experimental situation for CP violation in the B® —>■ chan¬ 

nel is not satisfactory, with measurements of the CP-violating asymmetries by the BaBar 
and Belle collaborations that are not in good agreement with one another. Measurements 
of CP-violation in the B Q S —>■ D~D J" channel might help to clarify the situation. In the 
case of the B ° —> D~Df mode, the LHCb collaboration has presented a first analysis 
of CP violation, with large experimental uncertainties. In the future, the experimen¬ 
tal errors can be significantly reduced. Using only information from branching ratios, 
we find the lower bound a > 0.052 for the penguin effects from the B® —>• and 

B® —> D~ Df branching ratios. Adding the current measurements of CP violation in 
B to the analysis, we obtain 

a = 0.35l°;“j , 9 = (215^)° , ^ = (60^)° 

from a y 2 fit to the data. These results indicate potentially sizeable penguin effects, 
although the large uncertainties do not allow us to draw further conclusions. 

Since the determination of from B® —>■ will not be competitive with the 

B—> J/'ipKg analysis and the control of penguins though B Q S —>• we advocate 

to use <pd as an input from the latter analysis for the determination of the penguin 
parameters from B® —> and relating them to their counterparts in B° -> D;Di 

with the help of the I/-spin symmetry of strong interactions. Following these lines, it 
will be possible to control the penguin effects in the determination of (j) s from the CP 
violation in the B ° —> D~Df channel. 

We find that the implementation of this strategy depends strongly on the values of 
the measured CP asymmetries, as we illustrated through a variety of future scenarios 
which are consistent with the current experimental situation, giving us a guideline for 
the LHCb upgrade era. We distinguish between two kinds of scenarios: in the first, the 
direct and mixing-induced CP asymmetries of the B® —>• channel are sufficient 

to determine a and 9 in a theoretically clean way, allowing us to determine the ratio 
|M'|/M| from the observable H, providing valuable insights into non-factorisable I/-spin- 
breaking effects. In the second — less favourable — class of scenarios, information 
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both from H and the CP asymmetries is needed to determine the penguin parameters. 
We have demonstrated that the resulting theoretical uncertainty for the penguin shift 
A (p s of the — > D~Df channel will be smaller than the experimental uncertainty for 
c t> e s s {B s —> DjD+) in both classes of scenarios. Similar analyses can be performed for 
B d — > D* d ~D d + and B ° — >• D*~D* + decays, where time-dependent measurements of the 
angular distribution of the decay products of the two vector mesons are required 77,78 


Analyses of B —> DD decays in the era of Belle II and the LHCb upgrade will offer 
interesting new insights both into the physics of strong interactions and into CP violation. 
We look forward to confronting the strategies discussed in this paper with future data! 
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Appendix: Notation 

In this Appendix, we give an overview of the notation, parameters and observables used 
in our analysis of the B —> DD system. 


Var. 

Eq. Amplitude ratio 

Description 

a 

X 

r P 

VA 

rpA 
r A 
'PA 
\„ |2 

( 

( 

I 

(43 

(42 

(65 

L02 

L85 

(41 

) 

) 

) 

) 

) 

) 

(P(“d + PA^)/(T+E+P^+PA^) 
(E + PAW)(T + PW) 

p(ct) jrp 

A/pH) 

PAH) j pH) 

(1 + r/ i)/(l + r PA ) 

~ T(B -A DD')/dT(B -> D^)/dq 2 

Penguin contribution wrt. total ampl. 
Exchange and penguin annihil. contr. 
Penguin contribution wrt. tree 
Annihilation contr. (in charged P’s) 
Penguin-annihilation contribution 
Comparison between A and PA 
Non-factorisable effects 

P 

£,<5 

( 

( 

(86 

L15 

L79 

) 

) 

) 


•SG^-breaking in T+P contributions 
S'G(3)-breaking in E+PA contributions 
SU (3)-breaking in a and 9 


Table 9: Overview of the various amplitude ratios and theoretical quantities. 
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Decay 

Amplitude 

Topologies 

T P E PA A 

Variables 

B° d D-D} 

A 

X X X X 

a, x, rp 

B° s -)• D~D+ 

A' 

X X X X 

a', xr' p 

B° s -A D~Df 

A 

X X 

a, x 

B° d -► D~ d Di 

A! 

X X 

a 1 , x 1 

B° d -> D-Df 

Aepa 

X X 

X 

B°, -► D d D+ 

A' 

EPA 

X X 

x' 


B+ - 

+ D ° D i 


X X 

X 

Q'ci'r A 

B+ - 

+ D°Df 

A' c 

X X 

X 

a'c,r' A 


Tabic 10: Overview of the variables relevant to the various decays. 


Decay Ratio 

Observable 

Eq. 

Value 

Used for 

B -A- DD'/B -A Dtv 

Rd 

( 

33|141) 


aiVF 

B ° d ^ D-D+/B" -A D-Df 
B° s -> D-Df/B° d -A D d Df 
B+ -A D°Df/B+ -A D°Df 

H 

H 

Be 

( 

( 

( 

25) 

129 

131 

) 

) 

( 

( 

( 

162 

163 

164 

) 

) 

) 

a 

a 

a c 

B°, -> D;Df/B« D d Di 
B° d -> O d D+/B°, DJD+ 


( 

( 

95) 

105 

) 

( 

( 

94) 

104 

) 

x! 

X 

B° s DfDf/B° d -A D~ d Df 
B° d -A D~D+/B° S -A D-Df 


( 

( 

117 

121 

) 

) 

( 

( 

116 

120 

) 

) 

X 1 

X 


Table 11: List of decay ratios and the corresponding observables and variables. 
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